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ABSTRACT 

AS AS, MACHO, OGLE and SAAO JHKL light curves of 13 stars, that have at some time 
been classified as D-type symbiotics, are analysed. Most of the near-IR light-curves that have been 
monitored over many years show long-term changes due to variable dust obscuration, in addition to 
the stellar pulsation. The distances to these objects are derived from the period-luminosity relation 
and estimates of the mass-loss rates made from the Kq — [12] colour. 

We reclassify AS 245 as an S-type symbiotic, with a semi-regular cool component with a pulsa- 
tion period of about one year. The periods of the large amplitude pulsations of SS73 38 (463 days), 
AS 210 (423 days) and H2-38 (395 days) are estimated for the first time, confirming that they are 
symbiotic Miras. 

A comparison of the symbiotic Miras with normal Miras of similar pulsation period shows that 
the symbiotic stars have on average higher values of A'o — [12] . This may indicate that they have 
higher mass-loss rates, or more likely that the dust which is being lost by the Mira is trapped within 
the binary system. 

Key words: binaries: symbiotic - Stars: individual: o Cet, RX Pup, V366 Car, BI Cru, SS73 38, 
V347 Nor, AS 210, AS 245, H 2-38, RR Tel, R Aqr, StHA 55, V335 Vul - sur-veys 



1. Introduction 

Symbiotic stars are long-period interacting binary systems, in which an evolved 
red giant transfers material onto its much hotter companion, which in most systems 
is a white dwarf. Based on their near-IR characteristics, symbiotic stars divide into 
two main classes (Allen 1982) depending whether the colours are stellar (S-type) or 
indicate a thick dust shell (D-type). The majority (~80%) of catalogued systems 
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are S-type and have near-IR colours consistent with cool stellar photosphere tem- 
peratures of ~ 3500 - 4000 K. Most of them have orbital periods ~ 500 - 1000 
days (e.g. Mikotajewska 2008). The near-IR colours of D-type systems indicate 
the presence of a dust shell which obscures the star and re-emits at longer wave- 
lengths. IR photometric monitoring has shown that these D-type systems have large 
amplitude variations and that they contain Mira variables with pulsation periods in 
the range 300-600 days; they are often called symbiotic Miras (Whitelock 1987). 
Since they must accommodate the Mira with its dust shell, these D-type systems 
should have much longer orbital periods than the S-types, a few tens of years and 
more. The latest review of symbiotic Miras and a comparison with normal Miras 
can be found in Whitelock (2003). 

Light curves of symbiotic stars reflect the very complex behaviour of these 
systems. They show high and low activity stages, flickering, nova-like outbursts 
originating from the hot component (S & D types), eclipses, ellipsoidal variability 
connected with orbital motion (S-type), radial pulsations (all D-type and some S- 
type) and semi -regular- variation (S-type) of the cool component, long-term dust 
obscuration (mostly D-type) and other types of variability (Mikolajewska 2001). 

In this paper we analyse the light curves of 13 objects that have at some time 
been classified as D-type symbiotics. The light-curves were provided by massive 
photometry surveys such as AS AS, OGLE, MACHO and near-IR monitoring at 
SAAO. 

2. Data 

Belczynski et al. (2000) listed coordinates for symbiotic stars, but many of 
these are not sufficiently accurate to identify the symbiotics unambiguously. So we 
first identified the 2MASS counterparts using the existing finding charts and the 
Aladin Java graphics interface running at the CDS in Strasbourg. This works well 
because symbiotic stars, which have the near-IR colours of late-type giants, are 
intrinsically bright in JHK. The 2MASS coordinates were then used to identify 
symbiotic stars in the OGLE, MACHO and ASAS databases. 

For o Cet, RX Pup, V366 Car, BI Cru, SS73 38, AS 210, RR Tel, R Aqr, 
StHA 55, and V335 Vul, light-curves were taken from the ASAS database (Poj- 
mariski 2002)0. These comprise V -band photometry obtained between November 
2000 and February 2009. The ASAS V light-curves are illustrated in Fig. 1, while 
Table 1 contains the ASAS names, the mean mags {V) and the full range of the 
variations found in V (AV). Although the ASAS database also contains SS73 38, 
there are only a few observations of it. 

The OGLE-II database (Udalski et al. 1997]! includes a light-curve for H2-38. 
This comprises /-band photometry obtained between 1997 and 2000. The OGLE 

'official home page of ASAS project: %ttp://www. astrouw.edu.pl/asas/\ 
^official home page of OGLE project: http://ogle. astrouw. edu.pl/\ 
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/ light-curve of H2-38 is shown in the second panel of Fig. 2, while Table 1 gives 
the OGLE name, the mean mag (7) and the full range of the variations found in / 
(A/). 

The MACHO database (Alcock et al. 1992)1 contains observations for H2-38, 
obtained between 1993 and 1999. The photometry was made through non-standard 
blue (Sm) and red (/?m) filters. The MACHO light-curve is shown in Fig. 2, while 
Table 1 gives the MACHO name, the mean mags (Bm» ^m) and the full range of 
the variations found in each band (ABm. A7?m)- 

We also analyse near-IR JHKL photometry for o Cet, RX Pup, V366 Car, 
BI Cru, V347 Nor, SS73 38, AS 210, AS 245, RR Tel and R Aqr. This was ob- 
tained with the 0.75-m, 1.0-m and 1.9-m telescopes at SAAO and is on the system 
described by Carter (1990). The photometry is illustrated in Figs. 3-5 and mean 
magnitudes, the full range of variations and the amplitudes of pulsation estimated 
by fitting sinusoids to detrended light curves (see Section 3.1..1 for period deriva- 
tion) are listed in Table 2. Some of the early data were published and discussed 
by Feast et al. 1983a, 1983b; Whitelock et al. 1983; Mikolajewska et al. 1999; 
Kotnik-Karuza et al. 2006, Santander-Garcia et al. 2007. The near-IR magnitudes 
of these objects are listed in appendix table and available in electronic form at Acta 
Astronomica Archive 0. 

Visual hght-curves for o Cet, RX Pup, BI Cru, AS 210, RR Tel and R Aqr 
were extracted from the database of the American Association of Variable Star 
Observer^ and used for comparison with the near-IR data. 

Most of objects discussed here are definitely symbiotic Miras, with four excep- 
tions; the symbiotic nature of o Cet, StHA 55 and V335 Vul is not certain. There 
have been suggestions (Jura & Helfand 1984, Kastner & Soker 2004, Ireland et 
al. 2007) that Mira B may be a low mass main-sequence (MS) star. StHA 55 and 
V335 Vul are only suspected of being symbiotic (Belczyriski et al. 2000) and may 
well be normal C-Miras (see also Munari et al. 2008). AS 245 is re-classifed here 
(section 3.1. .2) as S-type symbiotic star. These objects are indicated in figures by 
different symbols. 

3. Analysis and Results 

3. 1. Variability 
3.1..1 Period analysis 

All light-curves were analysed using the program perioeI ver. 5.0, based on the 
modified Lomb-Scargle method (Press & Rybicki, 1989). Long-term trends were 

^official home page of MACHO project: \ht tp://wwwmacho. mcmaster. ca/\ 
^official home page of Acta Astronomica: http://acta. ast rouw. edu.pl/\ 
^official home page of AAVSO: http://www. aavso. org/^ 
^the source program is available on http://www. starlink. rl.ac. uk/\ 



Table 1: D-type symbiotic stars in optical massive surveys: designation, average magnitudes and full range of the variations (A) 
during the observation run. The first column lists the identification number of the symbiotic star from Belczynski et al. (2000). 



No. 


Name 


Name in Survey 


V 


^v 


/ 


A/ Sm 
(mag) 






010 


o Cet 


ASAS 021920-0258.6 


6.75 


6.28 










026 


RXPup 


ASAS 081412-4142.5 


12.27 


1.53 










030 


V366 Car 


ASAS 095443-2745.5 


13.78 


1.95 










034 


BI Cm 


ASAS 122326-6238.3 


11.44 


2.90 










039 


SS73 38 


ASAS 125126-6460.0 














069 


AS 210 


ASAS 165120-2600.5 


12.91 


2.04 










120 


H2-38 


OGLE 180601.19-281704.2 
MACHO 105.21287.23 






12.21 


1.73 

-9.61 


>0.33 


-10.78 0.96 


175 


RRTel 


ASAS 200419-5543.6 


11.77 


0.42 










188 


R Aqr 


ASAS 234349-1517.1 


8.40 


5.86 










s03 


StHA 55 


ASAS 054642-^0643.8 


13.01 


2.90 










s26 


V335 Vul 


ASAS 192314-^2427.7 


11.80 


2.50 











Table 2: Near-IR SAAO photometry of the symbiotic stars: average magnitudes in the JHKL— hands, full range of the variation (A) 
of the JHKL magnitude and the amplitudes of pulsation estimated by fitting sinusoid with peak-to-peak amplitudes of ApJ etc to the 
observations. The first column lists the identification number of the symbiotic stai^ from Belczyriski et al. (2000). 































No. 


Name 


J 


H 


K 


L 


Ay 


A// 


aa: 


AL 


A/>y 


ApH 


hpK 


ApL 
















(mag) 












010 


o Get 


-1.23 


-2.10 


-2.55 


-3.12 


1.38 


1.40 


1.24 


1.39 


1.00 


1.01 


0.84 


0.68 


026 


RXPup 


5.67 


4.18 


2.98 


2.28 


3.72 


3.01 


2.30 


1.59 


1.08 


1.00 


0.81 


0.64 


034 


BICru 


7.62 


6.14 


4.84 


3.30 


2.07 


1.40 


1.11 


0.93 


0.69 


0.46 


0.31 


0.19 


030 


V366 Car 


7.19 


5.77 


4.78 


3.52 


3.05 


2.58 


1.93 


1.28 


0.53 


0.59 


0.49 


0.39 


039 


SS73 38 


9.41 


7.47 


5.96 


4.18 


3.18 


2.56 


1.79 


1.06 


0.53 


0.59 


0.49 


0.39 


060 


V347 Nor 


7.01 


5.54 


4.78 


4.02 


1.81 


1.32 


1.00 


0.94 


0.88 


0.70 


0.53 


0.44 


069 


AS 210 


9.50 


7.74 


6.28 


4.60 


2.67 


2.48 


1.84 


1.27 


1.06 


1.16 


1.05 


0.85 


105 


AS 245 


9.41 


8.05 


7.44 


6.92 


0.54 


0.57 


0.53 


0.44 


0.36 


0.42 


0.38 


0.30 


175 


RRTel 


6.57 


5.37 


4.43 


3.13 


3.47 


2.84 


2.05 


1.29 


0.71 


0.69 


0.60 


0.47 


188 


R Aqr 


0.80 


-0.27 


-0.89 


-1.72 


3.58 


3.07 


2.42 


1.79 


1.00 


0.91 


0.78 


0.72 



Vol. 58 



5 



Table 3: Pulsation periods derived from our analysis. 



No. 


Name 


Type 


ASAS 


OGLE 


MACHO 
(days) 


SAAO 


Other periods 


010 


o Cet 


O 


338±10 






332±3 


331.96' 


026 


RXPup 


O 


trend 






575±8 


578^ 


030 


V366 Car 





trend 






432±6 


433^ 


034 


BICru 





280±7 






274±2 


280'' 


039 


SS73 38 


c 








463±7 




060 


V347 Nor 











374±4 


373^ 


069 


AS 210 


c 


407±14 






423±7 




105 


AS 245 











366? 




120 


H2-38 







425±36 


395±16^ 






175 


RRTel 


o 


trend 






3 85 ±4 


385'' 


188 


R Aqr 


o 


395±13 






391 ±5 


386.96' 


s03 


StHA 55 


c 


372±15 








395'' 


s26 


V335 Vul 


c 


334±14 








347^ 



^ derived from ■ 

References: [1] Kholopov 1985; [2] Mikolajewska et al. 1999; [3] Feast et al. 1983b; [4] Whitelock 
et al. 1983; [5] Santander-Garcia et al. 2007; [6] Kotnik-Karuza et al. 2006; [7] Munari et al. 2008; 
[8] Sobotka et al. 2003. 



Table 4: New pulsation ephemerides. 



No. 


Name 


Ephemeris 




026 


RXPup 


Max{JHKL) 


= 2442238 + 575 x£ 


030 


V366 Car 


MaxiJHKL) 


= 2442413 +432 x£ 


034 


BlCru 


Max{JHKL) 


= 2443291 +274 xE 


039 


SS73 38 


Max{JHKL) 


= 2444983+463 xE 


069 


AS 210 


Max{JHKL) 


= 2446162 + 423 x£ 


120 


H2-38 


Max(.RM) = ' 


2449205 + 395 x£ 


175 


RRTel 


Max(JHKL) 


= 2442207 + 385 xfi 


s03 


StHA 55 


Max(V) =24522712 + 372 x£ 
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removed by subtracting a polynomial of appropriate order, and the resultant power 
spectra were compai^ed with the window specti'a. The periods were derived from 
the inverse of the maximum of the peak in the periodogram {P = f^ax)' whereas 
their accuracy was estimated by calculating the half-size of a single frequency bin 
(A/), centred on the peak {f^- is the center of the peak ) of the periodogram and 
then converted to period units {\P = f^^ - A/). The results of our period analysis 
are summarized in Table 3. Examples of our power spectra are shown in Fig. 6-7 
and Figs. 8 shows near-IR light curves folded with pulsation periods. 

The highest peak in a typical power spectrum corresponds to the pulsation pe- 
riod, while the other peaks represent annual aliases, second and third harmonics, 
long-term variation and some combination thereof. The only exceptions are the 
ASAS light-curve of V335 Vul, where the highest peak corresponds to the second 
harmonic, and the near-IR light-curves of RR Tel where it represents the annual 
alias. In both exceptional cases this is due to gaps in the light-curve of the object. 

There is practically no difference between the power spectra derived for the 
JHK or L observations. They all show the same position of peaks with a little 
difference in power. The latter is due to differences of amplitudes of the pulsations 
and long-term trends in the near-IR photometry. 

For four systems (AS 210, H2-38, SS73 38, and AS 245) the pulsation periods 
are derived for the first time. Pulsations are also detected in other systems with 
known periods. In particular, pulsations are always visible in the near-IR light- 
curves of all of the symbiotics we examined. However, the periods for o Cet and 
R Aqr (Table 3) ai^e not as accurate as those derived from visual data collected 
over a few centuries (Kholopov 1985), whereas for RX Pup, V366 Car, RR Tel 
and BI Cru our new estimates are better than published values (Mikolajewska et al. 
1999, Feast et al. 1983a, Whitelock et al. 1983) because they are based on more 
data. For some systems pulsations ai^e also visible in the optical light-curves, al- 
though the resulting periods are less accurate than those already known (Belczyriski 
et al. 2000, and references therein) because ASAS, MACHO and OGLE cover rel- 
atively short time periods (a few to several years). Finally, the optical light-curves 
of RX Pup, V336 Car and RR Tel do not show pulsations presumably because the 
Mira in these systems is obscured by optically thick dust and/or the optical light is 
dominated by emission from the hot component. 

The near-IR measurements for AS 245 suggest significant variability and a pe- 
riod of the order of one year, which we adopt for the following discussion. How- 
ever, with only 12 observations spread over almost two decades, the data are inad- 
equate for a proper analysis. 

Recently, Munari et al. (2008) estimated a pulsation period of 395 days for 
StHA 55 in the V -band. The value is very uncertain because they had only obser- 
vations covering only one pulsation cycle. The pulsation period of StHA 55 derived 
from ASAS data, which covers 6 pulsation cycles, is 372 days. 

The pulsations ephemerides for three systems for which periods are derived 
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here for the first time (AS 210, H2-38, SS73 38) and for StHA 55, RX Pup, 
V366 Car, RR Tel, and BI Cru are listed in Table 4. 

Fig. 9 presents the distribution of the pulsation periods for the symbiotic Mi- 
ras together with those for non-symbiotic Miras. For the symbiotic Miras we have 
combined the published periods (Belczyriski et al. 2000) with our new estimates 
(Table 3). The data for the non-symbiotic Miras are from Ohvier et al. (2001), and 
Whitelock et al. (1994, 2000, 2006). These non-symbiotic Miras are henceforth re- 
fen^ed to as 'normal Miras', but note that they will include widely separated binary 
systems; indeed o Ceti itself is included in both the symbiotic and normal group- 
ings. The normal Miras show distinct period distributions, with peaks at ~ 330 
and ~ 530 days for the O- and C-rich objects, respectively. The period distribu- 
tion for normal Miras is influenced by selection effect that are extremely difficult 
to quantify. We do know that there are 0-rich Miras with periods over 1000 days, 
the OH/IR stars, but these have progenitors of several solar masses and are quite 
rare. The symbiotic Miras have a mean period of about 400 days, and with one 
exception range from 280 (BI Cru) to 580 (RX Pup). The exception, V407 Cyg, 
has the period of 763 days (Kolotilov et al. 2003) and is thought to be hot bottom 
burning giant (Tatainikova et al. 2003) and is therefore probably more massive than 
its shorter period counterparts. There is no obvious difference between the periods 
of C-rich (5 objects) and 0-rich (19 objects) symbiotic Miras. All 0-rich symbi- 
otic Miras have periods longer than the median value for their normal counterparts. 
Whitelock (1987) suggested that this was essentially a selection effect - mass trans- 
fer in these systems is via the stellar wind of the Mira and long period Miras have 
in general higher mass-loss rates. Thus symbiotic activity will be seen from white 
dwarfs in binary systems with longer period Miras at much larger separations than 
from those in binaries involving shorter period Miras, i.e. we expect there to be 
unidentified white dwarfs in binary systems with short period Miras with very low 
levels of interaction. 

3.1..2 Long term trends 

The near-IR light-curves of all symbiotic Miras included in this study show, in 
addition to the Mira pulsations, significant long-term variations. Such trends are 
very common in symbiotic Miras, and are almost certainly caused by variable dust 
obscuration {e.g. Whitelock 1987; Mikotajewska et al. 1999). 

To get better insight into this phenomenon, we removed the Mira pulsations 
from the light-curves and examine the trends in the residuals. In the case of stars 
with periods longer than 400 days this is done by subtracting the best fitting sine- 
curve. For shorter period targets, o Cet, R Aqr, BI Cru and V347 Nor, the pulsation 
light-curves are asymmetrical and simply subtracting the sine-curve produced a lot 
of scatter in the residual. Therefore another method was used. We first generated 
an average pulsation curve and then subtracted it from the original light-curve. 
This method produced less scatter for the short period stai^s, while for the long 
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period pulsators both methods give similar results. In the case of BI Cru, data 
between JD2 445 400 and JD2 45 1 400 were used to prepare the average pulsation 
light-curves by excluding the dust obscuration events. In the case of R Aqr, data 
before JD2 444 000 were omitted from the average pulsation light-curves for the 
some reason. 

The light-curves, prior to removal of the pulsations, are plotted in the upper 
panels of Figs. 3 to 5. The secular trend is clearest at J and the behaviour oi J — K 
indicates that it is due to increasing reddening, as expected for increasing optical 
depth of the dust. This general trend is also present in the L light-curves. 

The near-IR colours of D-type systems indicate the presence of warm dust 
shells, bf Fig. 10 presents the J — K vs. K — L diagram for the symbiotic Mi- 
ras in this study together with those for normal Miras. The symbiotic Miras are 
shown in both their obscured and unobscured states. This demonstrates that most 
of the colours can be qualitatively reproduced with a shell of around 800K (see also 
Whitelock 1987). In at least 50 % of the studied systems the reddening toward the 
Mira star was larger than reddening toward the hot component {e.g. Mikolajewska 
1999). There appears to be very much less extinction towards the high excitation 
regions (emission lines and hot UV continuum) than towards the Mira suggesting 
that the hot component generally lies outside the dust cocoon associated with the 
Mira. This fact constrains the orbital separations in these systems to be > 10-15 
AU, and the orbital period to be > 20 yr, because the typical radius of a dust shell 
is > 5 RMira , and the Mira radius, RMira ~ 2-3 AU (Mikolajewska 1999). 

In the specific case of R Aqr the obscuration was explained in terms of orbitally 
related eclipses of the Mira by the accreting stream (Gromadzki & Mikolajewska 
2009, Willson, Gamavich & Mattel 1981). However, the dust obscuration phe- 
nomenon in symbiotic Miras cannot be, in general, orbitally related (Mikolajewska 
1999). First, these events seem to occur in well observed symbiotic Miras with too 
great a frequency, and in several systems with more than one event observed, they 
were separated by '-^ 2000 — 4000 days, much below the minimum orbital period 
expected for these systems. Secondly, while in RX Pup the infrared emission from 
the Mira shows modulation with a time scale of ~ 3000 days with two apparent 
minima around 1984 and 1990, at the same time the reddening towards the hot 
component was constant over at least 13 years (Mikolajewska et al. 1999). A ge- 
ometry of the binary system in which the cool component is eclipsed at least twice, 
while the hot one is not eclipsed at all is impossible. Moreover, the fact that no D- 
type system exhibits greater extinction towards the hot component than towards the 
Mira implies that only the Mira is affected. All these point to intrinsic variations in 
the Mira envelope as the main cause of the dust obscuration (see Sec. 3.2.2). 

For most of these systems the obscuration of light at J is not associated with a 
brightening at L , as might be anticipated if the dust forms in a uniform shell around 
the Mira (see Figs. 3 to 5). BI Cru is the one exception, but it has a spectrum quite 
unlike that of a normal 0-rich Mira in that it does not show the chai^acteristic H2O 
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absorption. What it does show is very strong 2.3 fim CO-band emission, and in 
that respect it has similarities to certain B[e] stai^s, such as Hen 3-1138 and Hen 3- 
1359, (Whitelock et al. 1983). It is therefore possible that we are observing a 
slightly different phenomenon in this case. We note that all parameters derived 
below (distance, mass-loss rate etc) for BI Cru through relations for normal Miras 
must be regarded as uncertain. 

C-rich Miras, including R For, also show obscuration events that cannot be 
reconciled with spherically symmetric dust ejection (Whitelock et al. 1997; Feast 
et al. . 2003; Whitelock et al. 2006). The alternative scenarios involve ejection 
around an equatorial disk or as puffs in random directions. In the cases of RR Tel 
and AS 210 spectro-polarimetry suggests a low inclination orbit (Schmid & Schild 
2002), but they show strong dust obscuration. This fact rather favours the ejection 
of puffs in random directions as was suggested for C-Miras (Feast et al. 2003; 
Whitelock et al. 2006) and is well known among R CrB stars. 

Secular trends are visible in the optical (ASAS) light-curves of RX Pup, V366 Car 
and RR Tel. These stars do not show pulsational variations in the visible and the 
trends observed are not well correlated with the long term changes in their near-IR 
light-curves. RR Tel and RX Pup underwent nova outbursts in the 1940s and 1970s, 
respectively, and the nearly linear trends in their optical light-curves are due to the 
fading of their hot components and nebular radiation. The cause of the optical trend 
in V366 Car is less clear, although the presence of a minimum with a time scale 
similar to the dust obscuration observed in its near-IR light-curve at earlier epochs 
may indicate some connection with a dust obscuration event. 

The visual light-curve of o Cet shows cycle-to-cycle changes of the visual am- 
plitude. Similar behaviour is seen in most Miras with sufficiently good light-curve 
coverage. Erratic variations ai^e also evident in the near-IR data, but the coverage of 
the light-curves is not as good. We note, however, that the highest visual maxima 
are correlated with the maxima in near-IR light-curve with (see left panels of Fig. 
3) which may suggest changes of the average luminosity of the Mira. 

AS 245 has colours typical of an oxygen-rich Mira (see Fig. 10). However, the 
amplitude of the possible pulsation (< 0.4 mag in all bands) is lower than in other 
symbiotic, and even non-symbiotic, Miras. Low amplitude and uncertain pulsation 
suggest that the cool component of AS 245 is an SRa variable rather than the Mira. 
We also point out similaiity between the cool giant of AS 245 and that in S-type 
symbiotic MWC 560 (V694 Mon). The latter also shows low-amplitude pulsations 
with a period of ~340'^ and low K — [12]. Gromadzki et al. (2007) suggested 
that the red component of MWC 560 is on thermally-pulsating AGB, although 
its pulsation characteristics may be influenced by its nearby companion possibly 
reducing the pulsation amplitude and the circumstellar dust. More observations ai^e 
necessary to settle the nature of the cool component of AS 245, but we have enough 
information here to reclassify it as an S-type; its position in Fig. 10 does not show 
signs of the dust associated with most symbiotic Miras. 
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3.2. Physical parameters 

The great advantage of long-period pulsating AGB stars is an opportunity to 
determine various physical parameters such as absolute magnitudes, distances and 
mass-loss rates using near-IR period-luminosity-colour relations derived by Feast, 
Whitelock and coworkers in a series of papers. In this section we apply these 
methods to estimate the physical parameters for symbiotic systems containing C- 
rich and 0-rich Miras. 

3.2..1 Extinction and distance 

The absolute K magnitudes of both O- and C-rich symbiotic Miras are estimated 
using the latest Whitelock et al. (2008) period-luminosity (PL) relationship: 

M^ = p[logP-2.38] + 5. (1) 

The slope p = — 3.51 ±0.20 was derived from large amplitude asymptotic gi- 
ant branch variable in the LMC. We use the zero-point 5 = —7.25 ±0.07 for O- 
rich Galactic Miras, estimated using the revised Hipparcos pai^allaxes together with 
published VLBI parallaxes for OH Masers and Globular Clusters. That value agrees 
with those estimated for 0-rich LMC Miras (5 = -7. 15 ±0.06) and C-rich Galac- 
tic (5 = -7.18 ±0.37) and LMC (5 = -7.24 ±0.07) Miras (assumes an LMC 
distance modulus of 18.39 ±0.05 mag; van Leeuwen et al. 2007). 

The interstellar extinction at V , Ay, is estimated using the Drimmel et al. 
(2003) 3-D Galactic extinction model, including the rescaling factors that correct 
the dust column density to account for small structures observed in the DIRBE data, 
but not included explicitly by the model. The extinction at JHK is then calculated 
from the relations given by Glass (1999) for photometry on the SAAO system. 
Obviously this extinction represents only the interstellar component and tells us 
nothing about any circumstellar reddening. 

We can compare the observed J — K with the intrinsic value, as discussed be- 
low, to get an estimate of the total, interstellar plus circumstellar extinction. The 
intrinsic {J — K)o for O Miras is estimated using the period-colour relation de- 
rived by Whitelock et al. (2000) for Miras in the solar neighbourhood observed by 
Hipparcos: 

(/_/i:)o = 0.71 log P- 0.39. (2) 

Intrinsic near-IR colours {J — K)o of C Miras are obtained from the period- 
colour relations for C Miras from Whitelock et al. (2006): 

(7-/:)o = 12.81 Hog P- 30.801. (3) 

Although the period colour relation for 0-Miras is quite well defined, at least 
at short periods, that for C-rich Miras is not; for example the best estimated rela- 
tion, {H — K)o vs. period, has a standard deviation of a = 0.48 mag. Therefore, 
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a reliable total extinction estimate is often impossible for C-rich objects and the 
uncertainty of the values derived for E{B — V) are huge. 

We assume wherever possible that the extinction derived from the period-colour 
relations is the total, circumstellar plus interstellar, value and that subtracting the 
interstellar extinction, evaluated as described above, gives the circumstellar value. 
The values of E{B — V) in the unobscured parts of light-curves, that are used for 
these estimates, are listed in Table 5. 

To estimate distances we derive the absolute K mag, Mk , from the period and 
take the observed K mag during intervals in which the objects did not show obvious 
dust obscuration, and average over the pulsation cycle. We use Ak derived from 
the colour excess to correct for the interstellar plus circumstellar extinction. The 
formal error of the distance derived by this approach is 12 - 20 %. The parameters 
derived as described above are listed in Table 4. 

3.2..2 Mass loss 

There is a good correlation between mass-loss rate and K — [12] colour for both 
C- and 0-rich Miras, provided we can assume that the shells are approximately 
spherically symmetric. This arises because to a first approximation the K and [12] 
mags are measures of the brightness of the star and of the shell, respectively. If, 
however, the shell is very asymmetric and our line of site to the star is obscured by 
thick dust (as it would be during an obscuration event) then we will underestimate 
the K brightness and overestimate the mass-loss rate. With this caveat in mind we 
can estimate the mass-loss rates for the symbiotic Miras from their K— [12] colours 
(Table 5). 

In the case of carbon Miras, Whitelock et al. (2006) estimated mass-loss rates 
using a modification of Jura's (1987) method, and fitted an analytical formula which 
we use here: 

logM = -7.668 + 0.7305(/:- [12]) -5.398 X 10-2^i5:-[12])2 (4) 
+1.343 xl0-^(ii:-[12])^ 

The mass-loss rates for the symbiotic 0-rich Miras are determined using the 
correlation with K — [12] colour derived for 58 high mass-loss 0-rich AGB stars in 
the South Galactic Cap (Fig. 21 of Whitelock et al. 1994). 

The [12] magnitudes are calculated from the IRAS 12-micron fluxes (Belczynski 
et al. 2000) using [12] = — 2.51ogFi2 + 3.62. The K magnitudes for the objects 
with SAAO near-lR light-curves are the values observed outside of obvious obscu- 
ration phases and averaged over the pulsation cycle. For HM Sge and V1016 Cyg 
we derived the average K from their published near-lR light-curves (Taranova & 
Shenavrin 2000). The use of a ^ magnitude obtained during a dust obscuration 
phase could result in an overestimate of the mass-loss rate by an order of magnitude 
(we assume, as suggested above, that the faint phases ai^e the result of asymmetric 
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obscuration and that the [12] mag remains approximately constant). For example 
the average K mag for V366 Car outside dust obscuration is 4.4 mag, from which 
we derive a mass-loss rate of 6.3 x lO^^Moyr^^; during dust obscuration K is 
about 6 mag and the derived mass-loss rate is 5 x lO^^MQyr^^ . The difference is 
not so large for objects with greater mass-loss rates, e.g., RX Pup outside of dust 
obscuration has M = 6.3 x lO^^M0yr^\ which increases to M = lO^^Moyr^^ 
during obscuration. 

For the remaining objects we use data from 2MASS (obtained between 1997 
and 2000) and Munari et al. (1992; obtained in 1990). Whereas the near-IR data 
in Munari et al. (1992) are from SAAO using the same photometric system as 
the light-curves presented here, the magnitudes from 2MASS have been trans- 
formed to the SAAO system using transformations from Carpenter (2001). A sig- 
nificant fraction of symbiotic Miras show dust obscuration events, so it is possible 
that some of these measurements were made during such events. Typically, for ob- 
jects with SAAO near-IR light-curves, the difference between the average K and 
transformed 2MASS Kg is 0.1-0.2 mag, and so we adopted the average of the Mu- 
nari et al. (1992) and transformed 2MASS values. The K magnitudes have been 
corrected for interstellar reddening (section 3.2..1). 

We should emphasize that mass-loss rates obtained in this way are not accu- 
rate. It is also important to appreciate that we do not yet understand the obscu- 
ration events or their link to the Mira mass-losing activity. If the symbiotic stars 
show obscuration events for the same reason that Feast (2003) and Whitelock et al. 
(2006) have suggested that the C-Miras do, then the mass-loss will be in discreet 
dust puffs in random directions like RCB stars (see also Whitelock 2003). The fact 
that symbiotic Miras are in binary systems will almost certainly effect what hap- 
pens to the dust once it leaves the immediate vicinity of the Mira in a non-random 
fashion, even if it does not affect the dust production itself. 

Fig. 11 presents the Kq— [12] colour distribution for symbiotic Miras and 
compai^es it to those for normal Miras. The normal 0-rich Miras show a distinct 
^0 — [12] colour distribution, which peaks at 1.76 mag, con^esponding to a mass- 
loss rate of ~ lO^^Mgjyr^^ , and has a tail that extends across the region occupied 
by the symbiotic Miras. As can be seen from Fig. 12, a plot of the — [12] colour 
vs. pulsation period, the tail is comprised of long period (mostly P > 400 days) 
objects. The C-rich Miras show a broad distribution from 1 to 9.5 mag which im- 
plies mass-loss rates ranging from ~ 10^^ to ~ 3 x lO^^Moyr^' . The symbiotic 
Miras show a distinct — [12] colour distribution with a peak at 4.25 mag, cor- 
responding to a mass-loss rate of 3.2 x lO^^MQyr^' . These figures show that 
the symbiotic Miras have on average lai^ger — [12] than the their normal coun- 
teipaits with similar- periods. Fig. 13 shows a plot of Kq — [12] vs. (/ — K)q for 
symbiotic and normal Miras. The scatter of symbiotic Miras is larger than that of 
the normal Miras and, although the distributions overlap, the symbiotic Miras have 
larger K — [12] than normal Miras of the same period or same J — K. If this is the 
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result of selection effects then we are only finding symbiotic Miras at the very end 
of their evolution when mass-loss is at its maximum. It seems more likely to be a 
consequence of the binary interaction. 

The widely accepted mass-loss scenario (e.g. Fleischer, Ganger & Sedlmayer 
1992) is that the Mira pulsation lifts matter above the atmosphere, but does not 
accelerate it to escape velocity. As the matter falls it encounters the next pulsation, 
which gives it another outward impulse. Matter is pushed by pulsations until its 
temperature drops enough for dust to condense (at a few stellar radii). Then, radi- 
ation pressure on the dust can efficiently accelerate the dust, and the dynamically 
coupled gas, to the escape velocity, — [12] is proportional to the amount of dust 
between us and the Mira, as described above. We should emphasize that in a binary 
system higher Kq — [12] could be the result of the lost mass not leaving the system 
rather than of intrinsically higher mass-loss rates. For example, if mass lost from 
the Mira in the binary was trapped within the system {e.g. near the L4 and L5 La- 
grange points), more cool dust in the system will mean more obscuration and more 
flux at 12 /urn and will result in higher values of Kq — [12] . 

An important question is about the influence of the secondary star on the mass- 
loss rate and its character. Podsiadlowski & Mohamed (2007) proposed a wind 
Roche lobe overflow (RLOF) model for o Get. In their model, a slow wind from 
Mira fills its Roche lobe and then the matter streams - via the LI Lagrangian point 
- onto an accretion disk around the companion. This models works most efficiently 
when the radius of dust shell is comparable with the Roche lobe radius (/?rl)- Or- 
bital separations for symbiotic Miras are essentially unknown, with the exception 
of R Aqr. The mean angular rotation rate from spectro-polarimetry suggest an aver- 
age period of around 150 yr (Schmid & Schild 2002). This estimation is uncertain 
because their observations cover about 10 yr (only about 5% of the implied orbital 
period). For typical component masses, Mhot ~ O.6M0 and MMira I.2M0, this 
period corresponds to an orbital separation of --^ 40 AU and 7?rl ~ 15 AU. Thus 
wind RLOF should occur in most symbiotic Miras. Podsiadlowski & Mohamed 
(2007) made SPH simulations, which show that in a case similar to symbiotic Mi- 
ras most of the material from the wind should be accreted or remain bound to the 
system, and only very small fraction is ejected to infinity, i.e., the accretion rate is 
~ 100 times higher than the Bondi-Hoyle value and the mass loss will tend to be 
confined to the orbital plane. 

4. Conclusions 

The study of symbiotic Miras is a difficult task because of the presence of an 
active accreting companion and the long time-scales of the variations (pulsation 
~l-2 years, dust obscuration --^10 years, orbital period ~100 yr). Fortunately 
data and relations estimated for normal Miras allowed us to derive some relevant 
physical parameters, such as distances, luminosities and order of mass-loss rates 



Table 5: The physical parameters of the symbiotic stars discussed here: distance (d), mean magnitude in the ^T— band outside of 
dust obscuration phases, absolute magnitude in K—band (Mk), average J — K colour (J — K), unredded J — K colour derived form 
colour-period relation {{J — K)o), interstellar- reddening estimated from a 3-D extinction model of Galaxy {Ej^_Y), total reddening 
estimated from colour excess (£'5_y), ^o-[12] colour and logarithm of mass-loss rate (logM). 



No. 


Name 


d 
(kpc) 


K 


Mk 


J-K 


(mag) 




Eb-v 


Kq-[\2\ 


logM 

(Meyr"') 


010 


o Cet 


0.11 


-2.55 


-7.75 


1.32 


1.40 


0.01 




3.04 


-6.2 


026 


RXPup 


1.6 


2.80 


-8.58 


2.36 


1.57 


0.51 


1.51 


4.90 


-5.2 


030 


V366 Car 


2.8 


4.40 


-8.15 


2.07 


1.48 


0.59 


1.13 


3.06 


-6.2 


034 


BICru 


2.3 


5.02 


-7.49 


2.58 


1.34 


0.91 


2.36 


4.25 


-5.5 


039 


SS73 38 


4.8 


5.34 


-8.27 


3.09 


3.40 


0.76 




3.68 


-5.6 


060 


V347 Nor 


2.9 


4.78 


-7.93 


2.23 


1.44 


0.49 


1.51 


2.83 


-6.3 


069 


AS 210 


5.6 


5.62 


-8.11 


2.72 


2.85 


0.27 




3.43 


-5.7 


105 


AS 245 


10.2 


7.44 


-7.89 


1.99 


1.43 


1.32 


1.03 


4.01 


-5.6 


120 


H2-38 


7.2 


6.60 


-8.01 


1.90 


1.45 


0.82 


0.85 


3.91 


-5.7 


175 


RRTel 


2.5 


4.18 


-7.97 


1.71 


1.45 


0.04 


0.51 


3.79 


-5.8 


188 


R Aqr 


0.24 


-1.06 


-7.99 


1.59 


1.45 


0.02 


0.23 


3.38 


-6.0 


s03 


StHA 55 


3.6 


5.27 


-7.92 


2.99 


2.13 


0.23 


1.64 


2.58 


-6.1 


s26 


V335 Vul 


3.7 


5.11 


-7.87 


1.85 


1.95 


0.50 




1.97 


-6.4 
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(provided certain assumptions are made). 

A comparison of the symbiotic Miras with normal Miras of similar pulsation 
period shows that the symbiotic stars have on average higher values of — [12] . 
This may indicate that they have higher mass-loss rates, or more likely that the dust 
which is being lost by the Miras is trapped within the binary system. Our obser- 
vations do not settle this issue. Undoubtedly, mass-loss rates and mass transfer in 
these systems are most interesting subjects, but more observational data, especially 
in UV, IR and radio bands, are needed to understand it. 
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Figure 1: ASAS light curves of type-D symbiotic stars. 
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Figure 2: MACHO (filled circle) and OGLE (open circle) light curve of H 2-38. 
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Figure 4: Near-IR light curves of V366 Car (left panels), SS73 38 (middle panels) and V347 Nor (right panels). 



o 



(N 



QC 

I 




8000 10000 12000 
JD-2440000 



RR Tel 1980 
. 2.5 



3.5 
5 



1990 



2000 



6 

A 

1 7 



^ 2 t 

3 
2.5 

_i 

3.5 
5 
6 

"= 7 



10 ■ 



12 







1 1 1 1 Illll 

*• V. If 

V ■ - 


Iiiiii 

•s; ■ 


1 1 1 hill 


1 : 

>*'V;v 


1 1 1 1'I'I 


Iiiiir 

; I.. 

• 


1 1 1 Illll 

• .'• • . ••• 


1 : 






' 1 1 1 Illll 

.*•..:. .jgLf ....... 

.1 . ■ TiiwiTisSi' L.^ • ' 


Iiiiir 








■ -tt': 



4000 6000 8000 10000 12000 
JD-2440000 



RAqr 


1 

A 

1 2 
3 
4 

-2 
-1 



1 

" 2 

3 

4 
5 



1980 1990 



2000 



> 


I II 


I 


1 1 1 [ 










» t 


















1 Ill 

^**^»*»#»« 


1 


i| 1 

»,.• 










i \ 

- VJ 


1 Ill 


1 

•*»'.•. 






1 Ill 

1 I* 


1 


:.,» 




















1 1 . ■ 


1 






4000 6000 8000 10000 12000 
JD-2440000 



Figure 5: Near-IR light curve of AS 210 (left panels), RR Tel (middle panels) and R Aqr (right panels). 
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Fie ure 6: Light curves of AS 210, H 2-38 and StHA 55 folded with pulsation 
periods (left panels) and related power spectra (right panels). Insights in top right 
corners of power spectra panels show power spectrum of windows. 
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Figure 7: Near-IR light curves of SS73 38 folded with pulsation periods (left pan- 
els) and related power spectra (right panels). Insights in top right comers of power 
spectra panels show power spectrum of windows. 
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Figure 8: Near-IR light curves of studied objects folded with pulsation periods. 
Ephemerides from Table 4 are used for most objects with the exception of o Cet 
and R Aqr. In case of these object ephemerides ai^e taken from Kholopov 1985. In 
panels of AS 245 dots represent SAAO data, triangle 2MASS, squares DENIS. 
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Figure 9: The pulsation period distribution for symbiotic Miras together with those 
for normal field Miras (Olivier et al. 2001, Whitelock et al. 1994, 2000, 2006). 
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Figure 10: The {J — K)o vs. {K — L)o colour-colour diagram. Colours are deredded 
using extinction from galactic 3-D model. Open dots represent oxygen rich sym- 
biotic Miras during dust obscuration, filled dots outside obscuration. Open squares 
represent carbon rich symbiotic Miras during dust obscuration, whilst filled outside 
obscuration. For comparison we also plot colours of oxygen rich Miras and car- 
bon rich Miras ( x ). For objects with uncertain nature different symbols were used: 
o Cet O), AS 245 (C). Straight line shows black body colours. Simple model of 
star with shell dust is plot as well. The temperature of star is 2750 K (black body), 
whilst dust shell temperature is 800 K. This model also includes line blanketing. 
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Figure 11: Kq— [12] distribution for symbiotic Miras together with those for nor- 
mal Miras. 



28 



A. A. 



CM 



7 
6 



4 



3 
2 
1 



« X 

x • • 

3 X X ^x X 

+ X • • X „ ^x 



+ + X 



x X ^ x^ 

X 

xx X ^ 

• X 



200 300 400 500 600 700 800 
P [days] 



Figure 12: Pulsation period vs. Kq — [12] colour for symbiotic Miras (•). For 
comparison we also plot colours of oxygen rich Miras (+) and carbon rich Miras 
(x). For objects with uncertain nature different symbols were used: o Cet (3), 
AS 245 (© ), StHA 55(A), and V335 Vul ( T ). The descrepant point at P=763 days 
is V407 Cyg, which is thought to be hot bottom burning giant and significantly 
more massive than the other symbiotic Miras (see text for detail). 
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Figure 13: {J — K)o vs. /Tq — [12] colour for symbiotic Miras (•). For comparison 
we also plot colours of oxygen rich Miras (+) and carbon rich Miras ( x ). For 
objects with uncertain nature different symbols were used: o Cet (3), AS 245 
(C), StHA55 (A), and V335 Vul(T). 



